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Abstract 

To investigate molecular epidemiology of dengue viruses (DENV) in southern China, a total of 14 dengue isolates 
were collected in southern China during each epidemic year between 1978 and 2006 and their full-length genome 
sequences were obtained by using RT-PCR method. The E gene sequences from additional 6 dengue fever patients 
in Guangzhou in 2006 were also obtained by using RT-PCR method. Combined with DENVs sequences published 
in GenBank, phylogenetic analysis and recombination analysis were performed. One hundred and twenty-five E 
gene sequences and 60 complete genome sequences published in the GenBank were also involved. Phylogenetic 
analysis showed that there was a wide genetic diversity of DENVs isolated in southern China. DENV-1 strains exist 
in almost all of the clades of genotype I and IV except the Asia 1 clade of genotype I; DENV-2 stains are grouped 
into four of the five genotypes except American genotype. DENV-4 strains are grouped into 2 genotypes (I and II). 
Phylogenetic analysis also showed that all DENV-4 isolates and two DENV-2 isolates were closely related to the 
prior isolates from neighboring Southeast Asia countries. The DENV-1 strain isolated during the 2006 epidemic is 
highly homologous to the strains isolated during the 2001 epidemic. 

Recombination analysis showed no inter-serotype recombination, but 22 intra-serotype recombination events were 
found across the 32 complete genomes of all Chinese isolates. The study suggested that dengue fever epidemic in 
Southern China over the past 30 years presented two important modes, 1) imported-cases-induced endemic 
prevalence; 2) endogenous epidemic outbreak with natural epidemic focus. Recombination may play an important 
role in dengue virus evolution and adaptation. 



Background 

Dengue fever (DF) and two more severe syndromes, 
dengue shock syndrome (DSS) and dengue hemorrhagic 
fever (DHF) are important mosquito-borne diseases in 
tropical and subtropical regions [1,2]. Since the first 
documented outbreak in 1779 in Jakarta, Indonesia, out- 
breaks have been documented in tropical and subtropi- 
cal regions. It has been the maximum public health 
burden in South-East Asia countries [3]. Dengue epide- 
miology changes fast among epidemic countries, and 
keeps on expanding to the non-epidemic area [1]. Since 
the pathogens were first discovered by Japanese scien- 
tists in 1943, dengue viruses (DENV) were isolated from 
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almost all South-East Asia countries including Thailand, 
Cambodia, Vietnam, Laos, Myanmar, Malaysia, Philip- 
pines, and Indonesia [3]. The dengue virus belongs to 
Flaviviridae family, and has four closely related but dif- 
ferent serotypes (DENV-1 through -4) in nature that are 
circulating or co-circulating [1,2]. Their genetic diversity 
remains wide, due to 1) absence of a proof-reading 
capacity in RNA-dependent RNA polymerases [4], 2) 
emergence of different lineages or clades during epi- 
demic [5,6], 3) increasing natural recombination [7-10], 
and 4) co-circulation of more than one serotypes in a 
locality [11-14]. 

Since imported DF epidemic had been reported in 
Hankou, Hangzhou, Shanghai and Guangzhou in 1920s 
and 1940s [15], there was no DF case reported in China 
till the outbreak occurred in 1978 in Foshan, Guang- 
dong Province, with DENV-4 infection [16]. And then 
DF was prevalent in Guangdong, Guangxi and Hainan 
province. A DENV-1 epidemic occurred in 1979 in 



O© 201 1 Wu et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons 
BiolVlGCl C^ntrBl Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in 
any medium, provided the original work is properly cited. 



Wu ef al. Virology Journal 201 1, 8:322 
http://www.virologyj.eom/content/8/1/322 



Page 2 of 9 



Zhongshan, Guangdong Province. This serotype of virus 
continued causing outbreaks over 2-3 years intervals. It 
become the dominating serotype and caused the latest 
outbreak in 2006 [17,18]. DENV-3 epidemic was only 
recorded once in 1980 at Zhan County, Hainan Island, 
where in 1985 DENV-2 caused an epidemic during 
which the first DHF case was reported [19]. After this 
DENV-2 epidemic, DENV-2 continued to be transmitted 
into Guangdong, Guangxi and Hainan until 2001, 
including three outbreaks in Foshan in 1993, 1998 and 
in Jiangmen in 2001 [20]. The second DENV-4 outbreak 
occurred in Guangzhou city in 1990 [20]. It was esti- 
mated more than 700,540 hospitalized cases with 513 
deaths from 1978 to 2007 [15]. Although the first isolate 
was sampled 30 years ago, dengue epidemic in southern 
China keeps increasing [15]. Lack of longitudinal 
research on dengue epidemics has hampered our under- 
standing of dengue molecular origin and evolution in 
China 

In this study, we determined the complete genome 
sequences of 14 dengue isolates collected in southern 
China during each epidemic year between 1978 and 
2006 and E gene sequences from six patients of Guang- 
dong, 2006. In combination with those published 
sequences in GenBank, we conducted an extensive 
molecular epidemiological analysis, aiming to determine 
where the DENV isolates in China originally came from, 
and what shaped their evolution. 

Materials and methods 

Ethics statement 

This research was approved by the Review Board of 
Guangzhou Medical Research Institute, the Review 
Board of Shenzhen Second people's Hospital, the Review 
Board of Beijing Institute of Genomics, the Review 
Board of Beijing Genomics Institute in Shenzhen and 
the Review Board of Southwest University. Written 
informed consent was obtained from the patient for 
publication of this case report and accompanying 
images. A copy of the written consent is available for 
review by the Editor-in-Chief of this journal. 

Sera 

Six dengue fever patients' sera were collected in 
Guangzhou city during the epidemic of DENV-1 in 
2006 by Guangzhou Medicine Institute. These sera were 
collected within 7 days after onset of symptoms and 
stored at -20°C. All sera were tested positive for DEN-1 
IgM by indirect immunofluorescence assay at the 
Guangzhou Medicine Institute. 

Viruses 

Fourteen DENV strains were obtained from patients 
sera in DF epidemic in China during 1978 to 2006 



Table 1 Description of DENV isolates sequenced in this 
study 



Serotype 


isolate 


Passage history 

(host, no. of passages) 6 


Locality 


■ ear 




GD03/91 


SMB, 3; C6/36, 1 


Guangzhou 


1991 




GD95/95 


SMB, 2; C6/36, 1 


Guangzhou 


1995 




GD01/97 


SMB, 1; C6/36, 1 


Chaozhou 


1997 


DENV-1 


GD99/99 


SMB, 1; C6/36, 1 


Chaozhou 


1999 




GD54/03 


SMB, 1; C6/36, 1 


Guangzhou 


2003 




GD66/03 


SMB, 1; C6/36, 1 


Guangzhou 


2003 




GD01/06 


SMB, 1 


Guangzhou 


2006 




GD02/06 


SMB, 1 


Guangzhou 


2006 




GD01/93 


SMB, 2; C6/36, 1 


Foshan 


1993 


DENV-2 


GD06/93 


SMB, 2; C6/36, 1 


Foshan 


1993 




GD08/98 


SMB, 1; C6/36, 1 


Foshan 


1998 




GD01/01 


SMB, 1; C6/36, 1 


Jiangmen 


2001 


DENV-4 


GD07/78 


SMB, 2 


Foshan 


1978 




GD09/90 


SMB, 2; C6/36, 1 


Guangzhou 


1990 



"All isolates were isolated from patients with primary infection manifesting 
classic DF symptoms at Guangdong Province, People's Republic of China. 
6 SMB, sucking mouse brain; C6/36, Aedes albopictus cell line. 



(Table 1). Sucking BALB/c mice and C6/36 cells were 
used to isolate viruses from clinical specimens. Sero- 
types of those isolates were confirmed by indirect 
immunofluorescence using anti-DENV monoclonal IgG 
(Guangzhou Medicine Institute, China). These isolates 
included eight DENV-1 isolates collected from Guangz- 
hou and Chaozhou City of Guangdong province in 
1991, 1995, 1997, 1999, 2003 and 2006 epidemic. Four 
DENV-2 isolates were collected from Fusan and Jiang- 
men City of Guangdong Province in 1993, 1998 and 
2001 epidemic. Two DENV-4 isolates were collected in 
Foshan City during 1978 epidemic and in Guangzhou 
city during 1990 epidemic respectively. 

Viral RNA extraction and RT-PCR methods 

Viral RNA was extracted using TRIzol reagent (Invitro- 
gen Corp., California, USA) according to the manufac- 
turer's instruction. The cDNA was synthesized by 
reverse transcription from 10 uL of RNA at 42°C for 50 
min in a 50 uL solution containing 50 mM Tris-HCl 
(pH 8.3), 75 mM KC1, 3 mM MgC12, 10 mM DTT, 100 
ng of the random hexamer primers, 200 U of Super- 
script II reverse transcriptase (Invitrogen), 25 U of RNa- 
sin (Invitrogen) and 0.5 mM dNTPs. The PCR was 
performed in a 25 ul mixture containing 2 ul of cDNA, 
10 mM Tris/HCl (pH 8.4), 50 mM KC1, 2.5 mM MgC12, 
100 uM each dNTP, 1 U Taq DNA Polymerase (Shang- 
hai Promage Corp., Shanghai, China), 0.25 uM each pri- 
mer (Additional file 1, Table S1-S3). The amplification 
reactions consisted of an initial denature step of 3 min 
at 95°C, followed by 35 cycles of 30 s at 94°C, 30 s at 
55°C, 60 s at 72°C. 
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Viral genome sequencing and E gene sequencing 

The complete genomes of the viruses were sequenced 
by a PCR-by-PCR strategy. At the same time, Envelope 
(E) gene from six patients' sera was also amplified. All 
primers were synthesized by BGI (Beijing Genomics 
Institute, Beijing, China) according to the reference 
sequence data of DENV-1 Singapore strain S275/90 
[GenBank: M87512], DENV-2 Jamaica strain 1409/88 
[GenBank: M20558], and DENV-4 Dominica strain 
814669/81 [GenBank: AF326573]. Amplified products 
were detected by agarose gel electrophoresis and 
sequenced using ABI 3730x1 DNA Analyzers (Applied 
Biosystems, Foster City, CA). 

Assembly of Genome Sequences and Sequence Analysis 

Genome assembly was performed independently by dis- 
tinct operators using Phred-Phrap-Consed [21]. The 
consensus sequence was yielded exactly with BLASTN 
and custom script which was written by Perl for all 
strains. A single contig was obtained for each of the 14 
isolates. The single contig was aligned to reference by 
BLASTN with the E-value of le-5. The computation of 
substitution tables were performed using a custom pro- 
gram written by Perl. 

Neighbor-joining trees were constructed using MEGA 
version 3.1 with the Kimura-2 parameter corrections of 
multiple substitutions [22]. Reliability of nodes was 
assessed by bootstrap resampling with 1,000 replicates. 
One hundred and twenty-five E gene sequences and 60 
complete genome sequences published in the GenBank 
were also involved. 

RDP3 was used for recombination analysis. We 
employed Bootscan, Chimera, GENECONV, MaxChi, 
and RDP methods incorporated in RDP3 Beta 27 pro- 
gram [23-27]. General recombination settings for all 
methods were as follows: Sequences were considered as 
linear, and the highest acceptable P-value was set to 
0.01, and event detected by two or more methods was 
taken into consideration. 

Results 

Nucleotide sequence accession numbers 

The determined DENV nucleotide sequences were 
deposited in GenBank database under the accession 
numbers FJ196841-FJ196860. 

Phylogenetic tree of DENV-1, -2, and -4 generated by E 
gene sequences 

To determine the DENV molecular epidemiology in 
southern China, 145 E gene sequences were used to 
infer a ML phylogenetic tree (Figure 1), including 1) 
20 E gene sequences determined in the study, 2) 
47 sequences with China localities available from 
GenBank, and 3) 78 representative sequences from 



GenBank referring to diverse genotype and geographic 
localities. 

The DENV-1 isolates from Guangdong belong to gen- 
otype I and IV (Figure IB). Within genotype I and IV, 
several sub-genotypes can be identified. Two isolates 
(CN.GD01/97 and CN.GD99/99) determined in the 
study clustered together as a single clade with a 96% 
bootstrapping support, which was named as a new Asia- 
4 clade. DENV-1 strains isolated in 2001 and 2006 
belong to Asia 3 clade in genotype I. However, the 
stains of 2006 could further be divided into three differ- 
ent sub-clades. Eleven strains of 2006 positioned at the 
end of the Asia 3 clade and the other four strains posi- 
tioned at two different sub-clades. Even the eight strains 
(two isolates and six patients' sera) collected in Guangz- 
hou city during the epidemic of DENV-1 in 2006 also 
clustered into 2 subclades in Asia 3 clade. All eight 
strains collected during 2002 and 2003 period at 
Guangzhou were identical and clustered closely as geno- 
type IV type. 

Chinese DENV-2 isolates were characterized into four 
genotypes except for American genotype (Figure 1C). 
Broadly, those collected in pre-1990s fell into Asia/ 
American and Asia II genotypes, while the post- 1990s 
belonged to Asia I and Cosmopolitan genotypes. Two 
isolates collected in 1998 in Foshan, which belonged to 
Asia I genotype, were identical to the strain from Thai- 
land in 1993. Three isolates of Jiangmen in 2001 epi- 
demic showed closely identical to the strain from 
Philippines in 2000 epidemic. 

Two DENV-4 epidemics occurred in China in the year 
of 1978 and 1990. Five isolates collected from the two 
epidemic years belonged to genotype I and genotype II 
(Figure ID). Two isolates collected during 1990 epi- 
demic belonged to genotype I and clustered closely with 
a strain isolated in 1984 epidemic of Philippines (PH. 
12123/84) [GenBank: U18435]. The three other DENV-4 
isolates from 1978 epidemic clustered together and fell 
into genotype II. They resided in phylogenetic tree 
much closer to the Indonesia isolate 30153 (ID. 30153/ 
73) [GenBank: U18428] collected in 1973 when there 
was a severe dengue epidemic. 

Phylogenetic tree of DENV-1, -2, and -4 generated by ORF 
sequences 

Based on the 74 complete coding regions (14 deter- 
mined in the study and 60 from GenBank), we con- 
structed maximum likelihood (ML) phylogenetic tree 
(Figure 2). 

In the study, the whole genome and the single E gene 
phylogenetic tree are in good coincidence with only lit- 
tle differences. Two DENV-1 isolates, CN.GD01/97 and 
CN.GD99/99, were assumed to be new members of a 
defined Asia-4 clade in E gene tree (Figure IB), while in 
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Figure 1 Neighbor-joining phylogenetic trees of 145 E gene sequences were conducted using Mega3.1 with the Kimura-2 parameter 
corrections of multiple substitutions. Bootstrap values (in percentage >65%) on each node were generated by using 1000 replications. The 
branch lengths are proportional to the number of nucleotide changes. All the Chinese isolates are in bold for clarity. A: phylogenetic tree of 
serotype DENV-1, DENV-2 and DENV-4 virus; B: phylogenetic tree of the DENV-1 isolates; C: phylogenetic tree of the DENV-2 isolates; D: 
phylogenetic tree of the DENV-4 isolates. 



the coding region phylogenetic tree; they clustered into 
the Asia-3 clade (Figure 2). 

Recombination detection 

RDP3 software was employed to examine the potential 
recombination events of all the Chinese isolates (43 
DENV-1,20 DENV-2,2 DENV-3,5 DENV-4) on E gene. 
Three recombinant events were detected on the E gene 
of DENV-1 with low P-value (P < 0.01) (Table 2), which 



was confirmed by phylogenetic construction (Figure 3). 
Neither DENV-2 nor DENV-4 has any evidence of such 
recombination across their counterparts. The potential 
recombination events were also analyzed across the 32 
complete ORF of all the Chinese isolates by RDP3 soft- 
ware. No inter-serotype recombination was detected, 
but 22 intra- serotype potential recombination events (17 
for DENV-1, 5 for DENV-2) were confirmed by more 
than two methods (Table 3). 
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Figure 2 Neighbor-joining phylogenetic trees of 74 complete coding regions were conducted using Mega3.1 with the Kimura-2 
parameter corrections of multiple substitutions. Bootstrap values (in percentage >65%) on each node were generated by using 1000 
replications. The branch lengths are proportional to the number of nucleotide changes. All the Chinese isolates are in bold for clarity. 



Discussion 

Since the first well documented DF outbreak in Foshan 
of Guangdong province in 1978, dengue has been 
reported periodically in southern China, with the latest 
epidemic at Guangdong province in 2006. In this study, 
14 Dengue isolates from 1978 to 2006 were sequenced 
and analyzed for molecular evolution. 



Dengue virus exists as four antigenically distinct 
viruses designated as serotypes (DENV-1 through 
DENV-4), belonging to genus Flavivirus of family Flavi- 
viridae. It has a positive-sense RNA genome that is 
translated as a single polyprotein and posttranslationally 
cleaved into three structural proteins and seven non- 
structural proteins [28]. The envelope protein (E) is 
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Table 2 Detected recombination events on the E gene 



region of Chinese DENV-1* 



Event 


Daughter 


Major Parent 


Minor Parent 


Region + 


1 


D1.GD23/95 


D1.GD14/97 


D1.GD95/95 


1-326 




(AY373427) 


(AY376737) 


(FJ 196844) 




2 


D1.GD01/97 


D1.NHD98039/98 


D1GZ02/03 


60-760 




(FJ 196847) 


(EF508198) 


(DQ2 11349) 






D1.GD99/99 






71-760 




(FJ 196848) 








3 


D1.GZ01/04 


D1.STD06098 


D1.71/02 


6-1299 




(EF032589) 


(EF508207) 


(EF025110) 





*Recombinant region and parental sequences were identified by five 
procedures (Bootscan, Chimarea, GENECONV, MaxChi, and RDP) integrated in 
RDP3 Package with a multiple comparison corrected P-value < 0.01. 
f Region coordinates are nucleotide positions of detected breakpoints in 
recombinant E gene sequence (1485bp). 



considered to be the immunodominant protein [29]. 
Dengue viruses also could be divided into different gen- 
otypes by the E gene [30,31]- So there is correlation 
between serotypes and genotypes of Dengue virus. Phy- 
logenetic trees of DENV-1, -2 and -4 generated by ORF 
and E gene sequences showed that there was a wide 
genetic diversity of DENVs isolated in southern China 
(Figure 1, Figure 2). DENV-1 strains exist in almost all 
of the clades of genotype I and IV but Asia 1 clade of 
genotype I (Figure IB, Figure 2). DENV-2 strains are 
grouped into four genotypes except American genotype 
(Figure 1C, Figure 2). DENV-4 strains are grouped into 
2 genotypes (I and II) (Figure ID, Figure 2). It also 
showed that different serotypes and genotypes epidemic 
prevalence exist in the same city and even during the 
same epidemic (Figure 1, Figure 2). It indicated that 
these dengue viruses maybe have different origination. 

Some of isolates, especially the first emerging of cer- 
tain dengue serotype strains in southern China were clo- 
sely related to those strains which were isolated in prior 
epidemics from neighboring Southeast Asia countries. It 
indicated that these dengue epidemics may be imported 
into China from the neighbor countries. For DENV-2, 
there were two outbreaks in Foshan in 1993 and 1998. 
Isolates from 1998 epidemic was closely related to the 
isolates from 1993 epidemic in Thailand (Figure 1C). The 
same situation was showed in DENV-4 isolates from 
1978 and 1990 epidemics. Isolates from 1978 epidemic 
was closely related to the isolates from 1973 epidemic in 
Indonesia (Figure ID). Isolates form 1990 epidemic was 
closely related to the isolates from 1984 epidemic in Phi- 
lippines (Figure ID). It indicated that these epidemics 
maybe imported into China from neighbor countries. 

Although some data showed the dengue fever were 
imported from neighbor countries, those epidemics with 
the same serotype of DENV during continues years 
showed that endemic infection of dengue circulating 
locally may be also the important cause of Dengue 



epidemic in southern China. In the study, two isolates 
(CN.GD01/97 and CN.GD99/99) clustered together as a 
new Asia-4 clade (Figure IB), suggesting that the iso- 
lates of 1999 maybe originated from the isolates of 
1997. Sequence analysis showed that the DENV-1 
strains isolated during the 2003 epidemic were closely 
related to the strains isolated during the 2002 epidemic 
of Guangzhou, suggesting that the 2002 strains had 
evolved in local and eventually cause the epidemic of 
2003 (Figure IB). Similar case was also showed in 
Guangzhou outbreak of DF in 2006. The DENV-1 strain 
isolated during the 2006 epidemic is closely related to 
the strains isolated during the 2001 epidemic (Figure 
IB). It suggested that DENV probably circulated in local 
and caused the epidemics of 2001 and 2006. Moreover, 
there were several clear divisions within the Chinese 
2006 isolates in the tip of the tree, even the eight strains 
collected in Guangzhou city by us during the epidemic of 
DENV-1 in 2006 also clustered into 2 subclades in Asia 3 
clade, showing viral evolution or increase of genetic 
complexity even during a single epidemic (Figure IB). 

Recombination plays a role in dengue virus evolution 
and adaptation [7]. We analyzed the potential recombi- 
nation events of all the Chinese isolates (43 DENV-1,20 
DENV-2,2 DENV-3,5 DENV-4) on E gene. Three 
recombinant events were detected on the E gene of 
DENV-1 with low P-value (P < 0.01) (Table 2). The 
incongruence of phylogenetic trees constructed sepa- 
rately by recombinant regions and non-recombinant 
regions confirmed the recombinant events (Figure 3). 
The recombination event 1 showed recombination 
between present GD14/97 isolate [GenBank: AY376737] 
as the major parent and the GD95/95 isolate [GenBank: 
FJ196846] as the minor parent, which led to a recombi- 
nant isolate of GD23/95 [GenBank: AY373427] (Figure 
3A). The recombination event 2 showed recombination 
between present NHD98039/98 isolate [GenBank: 
EF508198] as the major parent and the GZ02/03 isolate 
[GenBank: DQ211349] as the minor parent, which led 
to two recombinant isolates of GD01/97 [GenBank: 
FJ196847] and GD99/99 [GenBank: FJ196848] (Figure 
3B). The recombination event 3 showed recombination 
between present STD06098/06 isolate [GenBank: 
EF508206] as the major parent and the 71/02 isolate 
[GenBank: EF025110] as the minor parent, which led to 
a recombinant isolate of GZ01/04 [GenBank: EF032589] 
(Figure 3C). These studies of recombination provided an 
important context to understand the evolution of these 
viruses and sequence diversity generated in viruses in 
China. 

We also analyzed the potential recombination across 
the 32 complete ORF of all the Chinese isolates, and 
found no inter-serotype recombination occurred. The 
results are coincided with previous data [8,10,32]. 



Wu ef al. Virology Journal 201 1, 8:322 
http://www.virologyj.eom/content/8/1/322 



Page 7 of 9 



Non-i ecombination regions 



Recombination regions 



A 



i 



n 



m 



0005 



100 



,- CN.GZ01/02 DQ855296 



100 



100 
55 



-CN71AJ2EF025110 
„ 5 I — CN.GZD2T03 DQ211349 
CN.GZ01/95EF032590 
ACIMGD9995 FJ196846 
• CN.GD23f95AY373427 

CN.FJ231/04 DQ1 93572 

CN.NingBo/04 DQ83663 
CNGZ16/06EF113152 



CN.GD05/99AY376738 

OCN.GD14/97AY376737 



32M4S?bp 



001 



I 



100 



n 



CN.G7D1/02 DO 855296 
l_.CN.71/02EF025110 
tT"— CNGZ02/03 DQ211349 
A CN GD96T95 F J1 96846 



m 



100 



!^1CN GZ01/95 EF032590 
-CN.GZ16/06EF1 13152 

CN.FJ231/04 DQ1 93572 
CN.GD05/99 AY 376738 



51 
5-1 



H- CN.NingBo/04 DQ 83663 

♦ CNGD14/97AY376737 



t CN.GD23I95 AY37 3427 

l-326bp 



B 



0 005 



c 



m 



n 



m 



100 

— Lam 

i? A 



D1.GZ01/02DQ855296 
71/02 EF 025110 

D1.GZD2/03 DQ211349 
• D1.GD01/97 FJ1 96847 



a. 



55 



r— D1.FJ231/04DQ1 93572 

D1 NngBc/04 DQ836632 
D1.GZ17/06EF113153 
-♦D1.NHD98039/98EF508198 
62 j-« D1.GD99f99FJ196848 
_JD1.GD05/99AY376738 
D1.GD14/97 AY376737 



761-59bp 
761-^Obp 



ft 



CN.GZOUJ2DQ855296 
CN.7UD2EF025110 



7t'CN G2D2rt>3DQ21 1349 
.CN.GZ01/95EF032590 
H— CN.GD9SW5FJ196846 
ji,CN.FJ231rt)4DQ193572 
4 01 STD06098AK E F508206 
-•D1.GZ01/04EF032589 
:N.GD05/99AY376738 
CN.GD14/97 AY376737 

CNGZ16/06EF113152 
— CN.NingBortM DQ836632 




o o 1 



1300-5bp 



100 



r 



• D1.GD99/99FJ1 96848 



65 



D1.GZ01/02 DQ855296 
• D1.GD01/97FJ1 96847 
D1.71/02EF025110 
D1.GZ02/03 DQ211349 



m 



D1.FJ231/04DQ193572 

8£|— D1 NingBo/04 DQ836632 



■D1.GZ17/06EF113153 



♦ D1 NHD98039tf98EF508198 
— D1.GD05/99AY376738 
D1.GD14/97 AY376737 



0 005 



71-760bp 



100 



.-CN.GZ01/02DQ8S5296 



• D1.GZ01/04EF032S89 
J02EF025110 
GZ02/03 DQ211349 



n 



JQcN.71 
<J— CN.G 



CN G2D1/95 EF032590 



100 L CN.G095y95FJ1 96846 



39 



m 



100 



'5 



5m 



CN.FJ231/04 DQ193572 
♦ D1.STD06098/06EF508206 
CN.NingBo/04 DQ 836632 
CN. GZ1 6/06 EF 1131 52 
CN.GD05/99AY376738 
CN.GD14/97AY376737 



0 005 



6-1299bp 



Figure 3 Phylogenetic analysis based on non-recombination regions (Left) and recombination regions (Right) of event 1 (A), 2 (B), and 
3 (C) constructed by Mega 3.1 using Neighbor-joining methods. The event numbers are in accordance with table 2. Black circle, square, and 
triangle indicate the sequences of recombinant, major parent and minor parent, respectively. The trees are mid-rooted for clarity. 



However, we found 22 intra-serotype potential recombi- 
nation events (17 for DENV-1, 5 for DENV-2) (Table 3). 
Interestingly, for DENV-1, 13 of 17 potential recombina- 
tion events happened in the strains isolated in 1995, 
1997 and 1999 epidemics (Table 3), prior to which, in 



1995, the severest dengue outbreak ever reported, indi- 
cating a possibility of viruses escaping human immunity 
by recombination following epidemic. However, whether 
the presumption makes sense or just a coincidence 
remains an open question. 
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Table 3 Recombination events detected by RDP3 1 



Recombinant Event Number Breakpoint Positions in Recombinant Sequence(s) Minor Parental Sequence Major Parental Sequence 

Recombinant Sequence 





Begin 


End 








1 


677 


10104* 


D2.GD06/93 


D2.GD01/93 


D2.GD08/98 


2 


5514 


6548* 


D1.GD05/99 


D1.GD95/95 


D1.GD01/06 


3 


1* 


7487 


D1.GD01/97 


D1.GD14/97 


D1.GD95/95 








D1.GD99/99 


D1.ZJ01/04 


D1.GZ01/95 


4 


900 


1605 


D1.GD01/97 


D1.GZ71/02 


D1.GD14/97 








D1.GD99/99 


D1.GZ01/95 


D1.ZJ01/04 


5 


5763 


6401 


D1.GD99/99 


D1.GD66/03 


D1.GD01/06 


6 


1* 


466 


D1.GD05/99 


D1.GD23/95 


D1.GD99/99 


7 


5585 


6219 


D1.GD01/97 


D1.GD66/03 


D1.GD14/97 


8 


749 


1167 


D1.GD23/95 


D1.GD14/97 


D1.GD95/95 


9 


8240 


9343 


D2.GD08/98 


D2.GD01/93 


D2.GZ44/89 


10 


5457 


6201* 


D1.GZ71/02 


D1.GZ01/95 


D1.GD54/03 


11 


1606* 


2555 


D1.GD14/97 


D1.GD99/99 


D1.GD01/97 


12 


3711 


4913 


D2.GD01/01 


D2.GD06/93 


D2.ZS01/01 


13 


2588* 


2898 


D1.GD23/95 


D1.GD14/97 


D1.GZ01/95 


VI 


5171 


5506 


D1.GD23/95 


D1.GD14/97 


D1.GD95/95 


15 


7890 


8795 


D1.GD01/97 


D1.GD99/99 


Unknown 


16 


2708 


3121 


D1.GD01/06 


D1.GD54/03 


D1.ZJ01/04 


17 


4889 


5368 


D2.GD08/98 


D2.ZS01/01 


D2.GZ44/89 


18 


9645 


10146* 


D1.GD95/95 


D1.GD03/91 


D1.GZ01/95 


19 


7890* 


8385 


D1.GD99/99 


D1.GZ71/02 


Unknown 


20 


8820 


9106 


D1.GD99/99 


D1.GD23/95 


D1.GD05/99 


21 


2335 


2398* 


D1.GD95/95 


Unknown 


D1.GD23/95 


22 


1043* 


10104* 


D2.GD06/93 


D2.GD01/93 


D2.ZS01/01 


f = 22 final recombination events resulting from at least two method 
0.01 and one-by-one manual checking. Only the closest major and m 


i confirmation integrated in RDP3 package with a multiple comparison corrected P-value < 
nor parental sequences were shown in the table. Recombinants of DENV-1 were bolded for 



clarity. 

* = The actual breakpoint position is undetermined (it was most likely overprinted by a subsequent recombination event). 
Minor Parent = Parent contributing the smaller fraction of sequence. 
Major Parent = Parent contributing the larger fraction of sequence. 

Unknown = Only one parent and a recombinant need be in the alignment for a recombination event to be detectable. 



Conclusions 

Our study suggested that Dengue fever epidemic in 
Southern China over the past 30 years presents two 
important modes, 1) imported-cases-induced endemic 
prevalence; 2) endogenous epidemic outbreak with nat- 
ural epidemic focuses. Recombination analysis showed 
that recombination plays an important role in dengue 
virus evolution and adaptation. 
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